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We consider the longitudinal hydrodynamic evolution of the fireball cre- 
ated in a relativistic heavy-ion collision. Nonzero shear viscosity reduces 
the colling rate of the system and hinders the acceleration of the longitudi- 
nal flow. As a consequence, the initial energy density needed to reproduce 
the experimental data at RHIC energies is significantly reduced. At LHC 
energies, we expect that shear viscosity helps to conserve a Bjorken plateau 
in the rapidity distributions during the expansion. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

1. Introduction 

In heavy-ion collisions at relativistic energies a fireball of dense and 
hot matter is created. In later stages, the fireball expands and eventually 
hadronizes. In order to deduce the properties of that medium, including 
a possible phase transition to the quark-gluon plasma, a careful modelling 
of the evolution of the system is needed. With the assumption of local 
thermal equilibrium, the description of the dynamics of the system can be 
undertaken within the relativistic hydrodynamics [lj. Hydrodynamic cal- 
culations reproduce the spectra of particles in the transverse momentum, 
the collective flow, and the Hanbury-Brown-Twiss radii measured at Rel- 
ativistic Heavy Ion Collider (RHIC) [2]. Such calculations are essential in 
determining the correct equation of state of the matter under such extreme 
conditions, the initial conditions of the system and the freeze-out time. 
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Until recently, hydro dynamic calculations were performed assuming the 
applicability of the ideal fluid limit. It means that local equilibration pro- 
cesses are instantaneous, and no entropy is produced in the hydrodynamic 
stage. Relaxing the ideal fluid assumption amounts to the introduction of 
viscous effects in the evolution. For the fireball created at RHIC and Large 
Hadron Collider (LHC) energies, the most important dissipative effect in 
the hydrodynamic evolution is due to the shear viscosity [31 01 El [8] . A 
nonzero shear viscosity coefficient causes the saturation of the elliptic flow, 
a stronger transverse flow and can lead to a significant entropy production. 
Most of the calculations of the hydrodynamic evolution of the fireball with 
viscosity are done in the transverse directions only, with boost-invariance 
assumed in the beam (longitudinal) direction. However, shear viscosity is 
also important for the longitudinal expansion of the fireball [9]. In the fol- 
lowing, we present results of calculations in a 1+1 dimensional geometry 
of a non-boost-invariant expanding fluid with viscosity. The cooling rate 
and the acceleration of the longitudinal flow are reduced, and the entropy 
is produced in the expansion. 



2. Longitudinal expansion with viscosity 

Relativistic hydrodynamics with viscosity can be formulated consistently, 
without violating the causality [TO]. The hydrodynamic equations 

= (1) 

are modified; the energy-momentum tensor T^ u = TjfY , + tt^ u is composed 
of the energy-momentum tensor of an ideal fluid and a stress tensor tt„ v 
describing deviations from local equilibrium. Assuming the dependence 
of the densities only on the time t and the longitudinal coordinate z and 
restricting the flow velocity = (7, 0, 0, jv) only to the beam direction, 
one can write the hydrodynamic equations for a fluid with shear viscosity 
as [S] 

(e+p)DY = -ICp + UDY + )CU 
De = (e + p)KY - YIKY 

DU = (^ICY - n)/r w , (2) 

where D = u^d^ = cosh(Y - 8)d T + sinH ^ e) d e , K, = sinh(Y - 8)d T + 
coBhty-e) ^ = I to (f±§) is the space-time rapidity and Y = \ In (j^f) 
is the kinematic rapidity of a fluid element. Eqs. ([2]) involve four unknown 
function of 6 and of the proper time r = \/t 2 — z 2 : the energy density e, 
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the pressure p, the rapidity Y of the fluid element, and the shear correction 
II (in the 1+1 dimensional geometry the stress tensor tt^ v reduces to one 
independent scalar function II). One additional relation is given by the 
equation of state connecting e and p. We use a realistic equation of state by 
Chojnacki and Florkowski combining lattice QCD results and a hadron gas 
model |llj . The ratio of the shear viscosity coefficient r\ to the entropy of the 
fluid s is not known. We perform calculation for several values rj/s = 0.1-0.3 
of this parameter. For the relaxation time of dissipative corrections we take 
T n = 6r]/Ts [12] . The hydrodynamic equations ([2]) are solved in the r — 9 
plane starting from initial conditions at To = lfm/c 

e(r o ,0) = e o exp(-07(2<7 2 )) 

Y(t ,9) = 9 , (3) 

i.e. a Gaussian initial energy density profile and the Bjorken collective 
longitudinal flow of the fluid. 



3. RHIC results 

Eqs. ([2]) are solved numerically and the freeze-out hypersurface of con- 
stant temperature Tf = 165MeV is extracted. At the freeze-out hypersur- 
face, hadrons are emitted according to the Cooper- Frye formula [144 [3] 

dN visc dN d5N 

= 1 • (4) 

dy dy dy 

with the usual thermal emission contribution of the form 

dN S f 6max / , 
— = — { T (e) coshfy - 6) - t (6) sinhfy - 6 

dy 4vH J-e max V 

(2m£ + 2f + m 2 )i 

exp ^- mCOSh( ^-^ (g)) ) d9 , (5) 

where m is the meson mass, Yf(9) = Y(t(9), 9) is the fluid rapidity at the 
freeze-out hypersurface, and 

* = cosh(y-Y f (9)) ' (6) 
and an additional term due to the viscous corrections [1] 



d5N S 
dy 



/ ^ (r{9) cosh(y - 9) - r (9) sinh(y - 9)] 
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Fig. 1. Rapidity distributions for pions and kaons calculated with the ideal fluid 
hydrodynamics (dashed line) and using the viscous hydrodynamics with rj/s = 0.2 
(solid line) [5]. The dotted line denotes the results of the viscous hydrodynamic 
evolution, but neglecting the viscous corrections to the particle emission at freeze- 
out (Eq. El). Data are from the BRAHMS Collaboration [15] . 



At the freeze-out temperature, most of the pions and kaons come from 
secondary decays of heavier resonances. This effect is taken into account by 
a factor, equal to the ratio of all mesons to primary mesons[15j, multiplying 
the calculated distributions. The parameters of the initial distribution a 
and eo are adjusted to reproduce the pion and kaon rapidity distributions 
measured by the BRAHMS Collaboration [13]. With the increase of the 
shear viscosity of the fluid, two effects can be observed (Fig. |2J) in the 
retuned initial state of the fireball : 

• a reduced initial energy density of the fireball 



12£ 5 + 5£ 3 m 2 + 12£ 4 m + fm 3 - sinh(y - Y f {9)) 
(24£ 5 + 12£ 3 m 2 + 24£ 4 m + 4£ 2 m 3 + £m 4 ) 




(7) 
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Fig. 2. Initial energy density distributions for the ideal fluid hydrodynamic evo- 
lution with a realistic equation of state (dashed line) , for several viscous hydrody- 
namic evolutions (solid lines) , and for the ideal fluid with a relativistic gas equation 
of state (dashed-dotted line) [9]. All distributions after the hydrodynamic evolu- 
tion and freeze-out give pions distributions close to the BRAHMS measurements 

03- 

• and a small increase of the width of the initial distribution a with 
increasing viscosity. 

The two above mentioned effects are related to the change of the dy- 
namics induced by the shear viscosity. The reduction of the initial energy 
density at central rapidities is related to the reduced cooling rate in the 
viscous evolution. It originates from the smaller work of the fluid in the 
longitudinal viscous expansion, due to the change of the effective pressure 
from p to p — II (first equation in ([2D). The second effect is related to the 
acceleration of the longitudinal flow, as given by the second equation in ([5]). 
Gradients of the pressure in the longitudinal direction cause the acceleration 
of the flow, which becomes faster than the Bjorken one [16]. In the viscous 
evolution, the gradients of the pressure p are reduced by the gradients of the 
shear correction IT. As a result, the flow at the freeze-out is still Bjorken- 
like for the viscosity r//s = 0.2 (Fig. [3|), and is significantly slower than for 
the ideal fluid. Fast moving fluid elements emit hadrons far in the forward 
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Fig. 3. Difference between the flow rapidity of the fluid and the Bjorken flow taken 
at the freeze-out hypersurface, calculated for the evolution with the shear viscosity 
coefficient r\j s — 0.2 (solid line), for the ideal fluid with a realistic equation of 
state (dashed line) and for the ideal fluid with a relativistic gas equation of state 
(dashed-dotted line) [9]. 

and background rapidities, to counteract this effect a narrower initial energy 
density distribution in rapidity must be assumed for the ideal fluid. 

4. Predictions for LHC 

For the forthcoming LHC Pb-Pb heavy-ion experiments, we assume ar- 
bitrarily that the particle multiplicities at central rapidities would increase 
twice compared to Au-Au at the highest RHIC energies. This constraint 
serves to fix the initial energy density eo for the hydrodynamic evolution. 
At RHIC energies, experimental results show that a Bjorken scaling plateau 
at central rapidities, if existing at all, is very narrow [17\ [16| [9]. At LHC 
energies, we consider two different scenarios : a Gaussian initial energy 
density distribution in the space-time rapidity (similar as for RHIC), or a 
distribution with a plateau of width 2a p at central rapidities 

e(T O ,0) = e e W (-(d-a p ) 2 &(\9\ -a p )/2a 2 ) . (8) 
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Fig. 4. Initial energy density distributions for ideal fluid hydrodynamic evolutions 
(dashed and dashed-dotted lines) and for viscous hydrodynamic evolutions (solid 
and dotted lines). All distributions after the hydrodynamic evolution and freeze- 
out give the same pion density at central rapidity. 



The initial energy density distributions for the two scenarios are shown in 
Fig. HI both for the ideal fluid and for the viscous hydrodynamics with 
i]/s = 0.2. The initial energy density is smaller for nonzero viscosity. The 
reduction of the initial density accounts for the slower cooling and the en- 
tropy production in the later viscous evolution. 

The scenario with a Bjorken scaling plateau in the range of 8 units of 
central rapidities leads to different results than the one with the Gaussian 
initial density profile. Within the plateau region, the Bjorken scaling flow 
is stable, the pressure gradient in space-time rapidity is zero. During the 
evolution, the plateau region is becoming narrower, as the gradient of the 
pressure at the edges of the plateau starts to destroy the scaling form of 
the density and of the flow. If the shear viscosity is large, the rate at 
which the Bjorken scaling region is reduced is smaller than in the ideal 
fluid evolution. As a consequence, at the freeze-out, a substantial region 
with Bjorken scaling of the density and flow survives, if it is present in 
the initial state and if the shear viscosity is large (Fig. [5]). For the final 
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Fig. 5. Difference between the flow rapidity of the fluid and the Bjorken flow taken 
at t = lOfm/c, calculated for evolutions with shear viscosity coefficient rj/s = 0.2 
(solid and dotted lines) and for ideal fluid evolutions (dashed and dashed-dotted 
lines) for LHC initial conditions (Fig. [4|). 



meson distribution a plateau at central rapidities is possible only for the 
viscous evolution from an initial condition with a plateau in the initial 
energy density distribution (solid line in Fig. HJ). 

5. Summary 

We calculate the evolution of the matter created in relativistic heavy-ion 
collisions in the longitudinal direction both for the ideal and for the viscous 
fluid hydrodynamics Starting with Gaussian profiles of the energy den- 
sity in space-time rapidity and with a Bjorken scaling longitudinal flow, 
the hydrodynamic evolution reduces the energy density and accelerates the 
longitudinal flow. The rate of these processes is governed by gradients of 
the effective pressure p — IT, therefore nonzero shear viscosity corrections IT 
reduce the cooling rate and, at the same time, make the longitudinal flow to 
stay closer to the initial Bjorken scaling flow. At RHIC energies, comparison 
to meson rapidity distributions of the BRAHMS Collaboration [13] allows 
to constraint the parameters of the initial state. Depending on the value 
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of the shear viscosity coefficient one gets a reduction of the initial energy 
density by a factor 2-3 for rj/s = 0.1-0.2. At LHC energies, the same ef- 
fects take place for similar Gaussian initial conditions. Assuming an initial 
distribution with a plateau at central rapidities, where the Bjorken scal- 
ing solution applies, the dynamics is different. Within the plateau region 
the Bjorken flow is stable, both in the ideal and viscous hydrodynamics. 
Nonzero viscosity helps to preserve the Bjorken plateau in a wider region 
of rapidities trough the evolution; reduced gradients in the hydrodynamic 
equations make smaller the rate at which the Bjorken plateau diminishes. 
For r]/s = 0.2, the plateau region remaining till freeze-out is wide enough 
to survive the hadron emission process, and could be visible as a plateau in 
the final meson distribution in the kinematic rapidity. 



REFERENCES 

[1] P.F. Kolb and U.W. Heinz, Quark gluon plasma 3, edited by R. Hwa and 
X.N. Wang, World Scientific, Singapore, 2004, nucl-th/0305084. 

M. Chojnacki et al., (2007), arXiv:0712.0947 [nucl-th]. 

H. Song and U.W. Heinz, (2007), arXiv:0709.0742 [nucl-th]. 

D. Tcancy, Phys. Rev. C68 (2003) 034913, nucl-th/0301099. 

R. Baicr, P. Romatschke and U.A. Wiedemann, Nucl. Phys. A782 (2007) 313, 
nucl-th/0604006. 

R. Baier and P. Romatschke, Eur. Phys. J. C51 (2007) 677, nucl-th/0610108. 
A.K. Chaudhuri, Phys. Rev. C74 (2006) 044904, nucl-th/0604014. 
A. Muronga and D.H. Rischkc, (2004), nucl-th/0407114. 
P. Bozek, (2007), arXiv:0712.3498 [nucl-th]. 
W. Israel and J. Stewart, Annals Phys. 118 (1979) 341. 

M. Chojnacki and W. Florkowski, Acta Phys. Polon. B38 (2007) 3249, nucl- 
th/0702030. 

R. Venugopalan and M. Prakash, Nucl. Phys. A546 (1992) 718. 

BRAHMS, I.G. Bearden et al., Phys. Rev. Lett. 94 (2005) 162301, nucl- 
ex/0403050. 

F. Cooper and G. Frye, Phys. Rev. D10 (1974) 186. 

G. Torricri et al., Comput. Phys. Commun. 167 (2005) 229, nucl-th/0404083. 
L.M. Satarov et al., Phys. Rev. C75 (2007) 024903, hcp-ph/0606074. 
T. Hirano and K. Tsuda, Phys. Rev. C66 (2002) 054905, nucl-th/0205043. 



